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Abstract 

We derive a mass formula for n- dimensional unimodular lattices having 
any prescribed root system. We use Katsurada's formula for the Fourier 
coefficients of Siegel Eisenstein series to compute these masses for all root 
systems of even unimodular 32-dimensional lattices and odd unimodular 
lattices of dimension n < 30. In particular, we find the mass of even 
unimodular 32-dimensional lattices with no roots, and the mass of odd 
unimodular lattices with no roots in dimension n < 30, verifying Bacher 
and Venkov's enumerations in dimensions 27 and 28. We also compute 
better lower bounds on the number of inequivalent unimodular lattices 
in dimensions 26 to 30 than those afforded by the Minkowski- Siegel mass 
constants. 

1 Introduction 

First we review some definitions. For more information, see |9| and p2[ . 

An n-dimensional lattice A is the set Zvi + • • • + Zv n of all integer linear 
combinations of a basis {v±, . . . ,v n } for R n . We associate with A the Gram 
matrix A with i,j-th entry the inner product (vi,Vj). (The matrix A is a 
positive definite quadratic form, and much of what follows may be reformulated 
in the language of quadratic forms.) The determinant of A is defined to be 
det(j4), and A is called unimodular if its determinant is 1. We say that A is 
integral if (v, v') is an integer for all v,v' G A. Such a lattice is called even 
(or Type II) if (v, v) is always even, and odd (or Type I) otherwise. The dual 
of A is A' = {v E R™ : (v,x) € Z for all x £ A}, and the determinant of A is 
the order of the finite abelian group A'/A. Since we shall (with the exception 
of the occasional dual lattice) primarily be concerned with integral lattices, we 
generally omit the adjective integral in what follows. The norm of a vector 
v e A is defined to be (v, v), the square of its length. We say that a lattice is 
decomposable (or reducible) if it can be written as an orthogonal direct sum of 
two nonzero sublattices, and is indecomposable (or irreducible) otherwise. 

Each n-dimensional unimodular lattice A has a vector u such that (u, x) = 
(x, x) (mod 2) for all igA. Such a vector is called a parity vector ^ Preface 



to 3rd edition, p. xxxiv], or a characteristic vector Q, or a canonical element 
[ p3| . The set of parity vectors forms a coset u + 2 A in A/2 A, and each parity 
vector satisfies (u, u) = n (mod 8). 

Let N be an n-dimensional lattice with Gram matrix A and M be an Tri- 
dimensional lattice with Gram matrix B. We say that N is represented by M if 
there exists an to x n integral matrix X for which X t BX = A. We let r(M, N) 
denote the number of representations of N by M. When m — n, we say that 
M and N are (integrally) equivalent if there exists an integral matrix X with 
determinant ±1 for which X l BX = A. M and N are in the same genus if they 
are equivalent over the p-adic integers Z p for each prime p (including p — oo, for 
which Z p = Q). We define Aut(iV) to be the group of n x n integral matrices 
X for which X 1 AX = A. (Note that the definitions of det(iV), r(M,N) and 
equivalence are independent of the integral bases chosen for the lattices, and so 
is Aut(iV), up to isomorphism.) The theta series of a lattice A is defined by 
©a(<?) = E, 6 a 9 M = EZo KA, k)q k where q = e™. 

For our purposes, a root is a vector of norm 1 or 2, and the root system of A 
is the set of roots in A. The lattice generated by a root system is called a root 
lattice, and we define the rank of a root system to be the dimension of the cor- 
responding root lattice. Root lattices are completely classified: they are direct 
sums of the irreducible root lattices Z, A n (n > 1), D n (n > 4),Eq,Et, and E%. 
We use the same notation to refer to root systems, and for brevity we sometimes 
write the root system ri\A\ © • • • © njE 8 as A™ 1 ■ ■ ■ E^ 3 ■ 

In this paper we are concerned with the problem of classifying unimodular 
lattices, and also the subproblem of classifying unimodular lattices without roots 
(which correspond to denser sphere-packings), up to equivalence. 

The Minkowski-Siegel mass formula (see Q) gives the sum of the recipro- 
cals of the orders of the automorphism groups of all inequivalent lattices in a 
given genus. The mass constants can be used to verify that an enumeration of 
inequivalent lattices in a given genus is complete. They also give a lower bound 
for the number of inequivalent lattices in a given genus (sometimes called the 
class number). 

In dimensions divisible by 8 there are two genera of unimodular lattices, Type 
I and Type II (except in dimension 0, in which there is just one lattice of Type 
II) . In all other dimensions there is only one genus of unimodular lattices, Type 
I. The lower bounds provided by the mass constants show that the number of 
unimodular lattices increases super-exponentially as a function of the dimension. 
Unimodular lattices have been completely enumerated in dimensions n < 25, 
but number more than 900 million in dimension 30. Unimodular lattices without 
roots have been completely enumerated in dimensions n < 28, but number more 
than 8 x 10 20 in dimension 33. In both cases, somewhere between dimension 26 
and dimension 32 there is a transition from being completely classified to being 
numerous enough to make classification unappealing, so information about what 
is going on in these dimensions is of interest. 

Our approach is based on a suggestion by Borcherds that the Fourier coeffi- 
cients of Siegel Eisenstein series of degree 4fc could be used to derive something 
analogous to the Minkowski-Siegel mass formula, but which gives the mass of 
all the even unimodular lattices of dimension 8fc having any given root system. 

We wrote a computer program which uses Katsurada's formula (see JL6| ) 
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for the Fourier coefficients to calculate the masses for all possible root systems 
of even unimodular 32-dimensional lattices. From these masses, we used the 
methods of j^, Chapter 16] to compute the masses, for each root system, of 
unimodular lattices in dimensions n < 30. 

In particular, we have a mass formula for those lattices which have the 
empty root system (that is, which have no roots). This formula verifies the 
known masses in dimensions n < 28, and provides new lower bounds for the 
number of unimodular lattices without roots in dimensions 29 to 32. 

We also used our program to compute better lower bounds on the total 
number of odd unimodular lattices in dimensions 26 to 30, and even unimodular 
lattices in dimension 32, than those gotten from the Minkowski-Siegel mass 
constants. 

Our results may be viewed as a coarse classification of lattices, as inequivalent 
lattices may have the same root system. For dimensions n < 23, in which it 
happens to be the case that a unimodular lattice in a given genus is completely 
determined by its root system, our results coincide with the previously known 
enumerations. The same is true for the even unimodular 24-dimensional lattices 
(which are known as the Niemeier lattices). 



2 A mass formula for even unimodular lattices 
having any given root system 

Let SI be the set of inequivalent even unimodular lattices of dimension 8k. We 
define the mass m of SI by 

w =E|AuT(A)|- 
By the Minkowski-Siegel mass formula, for k > we have 

_ _ \B ik \ 4 -rr 1 \B 2 j\ 



n 



where Bi is the ith Bernoulli number. (See ji| and |jj for this, and for the mass 
formulae for other genera of lattices.) 

For an n-dimensional lattice N, we define a weighted average number of 
representations of iV by the lattices A £ SI by 



a(N) = 1 £ 

711 ' * 



r(A, N) 



m A^J Aut ( A )l' 

Let {Ri, ...,R S } be the set of all the root lattices of dimension n or less with 
no vectors of norm 1. Each Ri is the direct sum of lattices Aj(j > l),Dj(j > 
4:),Eq,Et, and Eg. Let J!,; be the set of lattices in SI having root system Ri, 
and let m(Ri) = J2Aen l^ ut (^)l _1 be the mass of those lattices in SI which 
have root system Ri, so that m = m(Ri) + • ■ ■ + m(R s ). 
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Proposition 1 Let U be the sx s matrix with i, j-th entry r(Rj , R{), let v be the 
vector (m(i?i), . . . , m(R s )) t ) and let w be the vector (a(i?i), . . . , a(i? s ))*. Then 
^Uv — w. Furthermore, U is invertible, so m(-Rj) = (mU~ 1 w)i gives the mass 
of the lattices in f2 having root system Ri . 

Proof. Observe that if R is any root lattice and if S is the root system of a 
lattice A then r(A, R) = r(S, R), so for j = 1, . . . , s we have 



i=l Aefii 1 

1 s 

m ^ — ^ 



m 

i=l 



Thus — Uv — w. We may assume the i?,'s are ordered so that their dimen- 
sions are non-decreasing, and so that among those with the same dimension 
the determinants are non-increasing. With this ordering, r(Rj,Rj) = when- 
ever i > j, so U is upper triangular. Since each diagonal element r(Ri,Ri) = 
|Aut(i?j)| is positive, we have det({7) ^ so v = mU~ 1 w. 



Remark 2 The values a(Rj) are the Fourier coefficients of Siegel Eisenstein 
series. Borcherds, Freitag and Weissauer Q used a relation similar to -^Uv = w 
to compute the coefficients of a cusp form from the known values of m(Ri) in 
dimension 24. We shall do the inverse in dimension 32: use the values of a(Rj) 
to derive the values m{Ri). We discuss how to compute a(Rj) in Section 0. 
The problem of computing r(Ri,Rj) is largely combinatorial; we discuss it in 
Section |[ 



3 Masses of 32-dimensional even unimodular lat- 
tices with any given root system 

Let mj/ (R) and m l n (R) denote the masses of the n-dimensional even and odd 
unimodular lattices having root system R, and let m n {R) denote their sum. 
Let w(R) denote the order of the Weyl group of R; w(R) is the product of the 
orders of the Weyl groups of the irreducible components of R, where w(A n ) = 
{n+\)\,w{D n ) = 2 n - 1 n\,w(E 6 ) = 2 7 3 4 5 • 7,w(E r ) = 2 10 3 4 5 ■ 7, and w(E s ) = 
2 14 3 5 5 2 7. It is sometimes more convenient to list the values of (R) ■ w(R) 
than it is to list the values of m^ 1 (R) alone, but the latter can easily be recovered 
from the former. 

We used a computer to calculate (R) for each root system R with 
rank(i?) < 32. The computation took about two weeks on a Sun Ultra 60, 
running a program written in Common Lisp and compiled with Franz Inc.'s Al- 
legro CL. (We discuss several issues related to the implementation in Sections 
to |.) 
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Theorem 3 The nonzero values of m^ 2 (R) ■ w(R) with rank(R) < 9 are as 
listed in Table |. (The list of all 13218 nonzero values of m^ 2 (R) ■ w(R) is 
available electronically at pl\j.) 

The root system R of a lattice A is called complete if the sublattice of A 
generated by R has finite index in A, or equivalently if rank(i?) = dim(A). The 
classification of even unimodular lattices with complete root systems is closely 
related to the classification of certain self-dual codes (see |34|). 

Corollary 4 (Kervaire |is|l ) There are 119 complete root systems R which 
occur as root systems of indecomposable 32 -dimensional even unimodular lat- 
tices. 

Proof. There are 143 root systems R for which rank(i?) = 32 and m^ 2 (R) > 0. 
Of these, the root system D\ & and the 23 root systems containing Eg corre- 
spond to decomposable lattices. The remaining 119 root systems correspond to 
indecomposable lattices. 

Kervaire |[8[ , extending the work of Koch and Venkov [^5|, has proven a 
stronger result: there are exactly 132 indecomposable even unimodular 32- 
dimcnsional lattices with complete root systems, with 119 different root systems 
occurring. 

Corollary 5 There are at least 10000000 32 -dimensional even unimodular lat- 
tices without roots. 

Proof. Each such lattice has at least two automorphisms, so the number of such 
lattices is at least 2 x m^ffl) = 1-096 x 10 7 . We will have more to say about 
lattices without roots is Section |H. 

Remark 6 If there is a single even unimodular 32-dimensional lattice with root 
system R then the corresponding value m 1 ^ (R) -w(R) is of the form 1 jq for some 
q <E N; the converse often holds (but not always: for R = A[A 2 A^Dj, m l ^ 2 (R) ■ 
w(R) splits as 1/4 = 1/12 + 1/6; see Q). We give one example below, and shall 
see more examples in Section ^ (Niemeier lattices) and Section || (unimodular 
lattices with no roots). The last column of Table [l] labels as "unique" those 
lattices whose uniqueness follows from [|J , and in some cases gives the name of 
the lower dimensional lattice which when glued to R produces the unique even 
unimodular 32-dimensional lattice with root system R. 

Example 7 In |l3| , Elkies and Gross construct a 26-dimensional even lattice Lq 
of determinant 3 with no roots, for which |Aut(L )| = 2 13 3 5 7 2 13 = 1268047872, 
and refer to a preprint fl2]| containing a proof of its uniqueness which uses 
Euclidean lattices. (Borcherds Chapter 5.7] had previously proved its ex- 
istence and uniqueness using Lorentzian lattices). By Q], even 26-dimensional 
lattices of determinant 3 with no roots are in one-to-one correspondence with 
even unimodular 32-dimensional lattices having root system Eq, where the or- 
der of the automorphism group of the 32-dimensional lattice is w{E§) times 
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Table 1: Masses of 32-dimensional even unimodular lattices 



dim root system R 



mg (R) ■ w(R) 



(as decimal) comments 








131003733 12820233266589 17 
238863431761920000 


5484461.50 


1 




1 1 1 E9<31(!01 UOI'IO 
1110iiDlDol»z4oo 

5677056 


19646832.00 


2 


i 2 

At 


72024731351193941 
1857945600 


38765792.00 


2 


A 2 


1327104974887 
2939328 


451499.44 


3 


t 3 

Ai 


6904800898075 
124416 


55497696.00 


3 


AiA 2 


QTT'Q^I 9^1 9 "37 

445440 


2195472.50 


3 


A 3 


ozy Iz rao 1 
74240 


4433.30 


4 




30223371257980501 
471859200 


64051672.00 


4 




1 osi(;7i m 
lytto ( hjimio 

3456 


5748582.50 


4 


) 2 


1772535692573 
42598400 


41610.38 


4 




z ± u 1 ( 

768 


27439.89 


4 


A 4 


8397751 
384000 


21.87 


4 


£>4 


35841940559 
157212057600 


0.23 


5 


A{ 


14457125482723 
230400 


62747940.00 


5 


A\A 2 


lUD.tfDoo4J.3l J 1 P>A 

995328 


10676264.00 


5 


A A 2 


673556587 
2560 


263108.06 


5 


A? A 3 


9nn'-i«f;«n'^7nQ 

Z4 JU. >.M).MJ..1 1 U!.' 

2211840 


90597.34 


5 


A2A3 


8085187 
5760 


1403.68 


5 




AfiQI 789^ 

269568 


174.05 


5 


A1-D4 


473 
240 


1.97 


5 


As 


73 
960 


0.08 


5 


As 


433 
3317760 


0.00 


6 


At 


6635520 


53604776.00 


6 


A A A 

A\A 2 


22484458507 
1440 


15614208.00 


6 


A 1 A 2 


16820220686833 
19169280 


877457.06 


6 


A 3 ,» 


25920 


210345.19 


6 


A| 


9785018477 
1866240 


5243.17 


6 


A1A2A3 


680633479 
61440 


11078.02 


6 


A\Ai 


11264777 
15360 


733.38 


6 


A 2 


323400013 
20447232 


15.82 


6 


A?D 4 


622763 
69120 


9.01 


6 


A 2 A 4 


6059 
512 


11.83 


6 


A2D4 


344077 
2446080 


0.14 


6 


A1A5 


123853 
159744 


0.78 


6 


A1D5 


1 

048 


0.00 


6 


A 6 


1 

4608 


0.00 



no roots 
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Tabic 1: Masses of 32-dimensional even unimodular lattices (cont.) 



dim root system R m 32 (R) ■ w(R) (as decimal) comments 



6 


D 6 


6 


Ea 


7 


A\ 


7 


A\A 2 


7 


A\A\ 


7 


A\A 3 


7 


AxA\ 


7 


A\A 2 A 3 


7 


A\A A 


7 


A\A 3 


7 


A\D A 


7 


A^A% 


7 


A 1 A 2 A i 


7 


A\A 5 


7 


A X A 2 D 4 


7 


A 3 A 4 


7 


A 2 Ar, 


7 


A 3 D 4 


7 


AjD 5 


7 


AxAe, 


7 


A 2 D 5 


7 


A 7 


8 


Af 


8 


A\A 2 


8 


A\A\ 


8 


A\A 3 


8 


A\A\ 


8 


A\A 2 A 3 


8 


A\ 


8 


A\A A 


8 


A X A\A 3 


8 


AfD 4 


8 


A\A\ 


8 


A\A 2 A 4 


8 


A 2 A\ 


8 


A\A 5 



18720000 
1 

1268047872 
23307759701 
576 

290429642677 

15360 
32712111919 

16128 
491983248817 
1290240 
641243179 



15360 
528426689 
11520 
828763 
384 
238819303 
552960 
2242333 

77760 
2750969 

17920 
57178267 
483840 
24071 
5760 
959 
640 
297043 
860160 
51383 
774144 

133 
30720 

11 
1152 
521 



188160 
1 

6912 
1 



1376256 
2005621383142854931 
73987522560 
2015000372681 



103680 
1327413084613 
368640 
980176289 
1728 
13761649541 

80640 
42164571593 
322560 
176324027322323 
180592312320 
9542951321 

1935360 
691208023 

161280 
371231029 



5160960 
6071138573 
7741440 
1655917 

2688 
1794539 
120960 
2245489 



0.00 
0.00 

40464860.00 
18908180.00 
2028280.70 
381311.40 
41747.60 
45870.37 
2158.24 
431.89 
28.84 
153.51 
118.18 
4.18 
1.50 
0.35 
0.07 
0.00 
0.01 
0.00 
0.00 
0.00 
27107562.00 
19434804.00 
3600838.50 
567231.70 
170655.37 
130718.54 
976.37 
4930.84 
4285.76 
71.93 
784.24 
616.04 
14.84 
15.66 



unique 
unique 
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Table 1: Masses of 32-dimensional even unimodular lattices (cont.) 



dim root system R m 32 (R) ■ w(R) (as decimal) comments 

8.49 
5.93 
4.15 
0.07 
0.86 
0.04 
0.05 
0.02 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 

0.00 unique 
0.00 odd Leech 



8 


A\A 2 D 4 


8 


A\A 4 


8 


AjA 3 A 4 


8 


A\D 4 


8 


AxA 2 A h 


8 


A{D b 


8 


A 1 A 3 D i 


8 


a\a 6 


8 


A\ 


8 


A z A h 


8 


A X A 2 D,, 


8 


A 2 A 6 


8 


A 4 D 4 


8 


Dl 


8 


A :i D 5 


8 


AiD 6 


8 


AU 7 


8 


A 8 


8 


A 2 E 6 


8 


D s 


8 


E$ 


9 


Al 


9 


A\A 2 


9 


A\A\ 


9 


A\A 3 


9 


A\A\ 


9 


AiA 2 A 3 


9 


AiAi 


9 


A\A 4 


9 


A\A\A 3 


9 


A\D 4 


9 


A\A\ 


9 


A\AiA A 


9 


A\A Z 


9 


AfA 5 


9 


AlAiDi 



6113 
720 
237 
40 
200995 
48384 
18917 
266112 
7739 
8960 
829 
20160 
353 
6720 
2287 
120960 
432673039 
170311680000 
1049 
483840 

43 
20160 
1 

3780 
1 



15360 
1867 



1937768448 
19 
1935360 

607 
85155840 
5 

1064448 
1 

185794560 
1 

277136640 
1 

1002795171840 
1 

8315553613086720000 
4668288705497 

290304 
5526347655971 



322560 
740728763 
144 

1568559349553 

2211840 
1368567316381 
2903040 
1271443287 
4480 
23738023 
2560 
2142275851 

233280 
927969995 
43008 
158461 



1080 
36859871 

13440 
237492443 
107520 
3078967 
17280 
457601 



10080 
21454849 






.00 


16080690. 


00 


17132774. 


00 


5143950. 


00 


709164. 


94 


471425. 


60 


283804. 


30 


9272. 


67 


9183. 


.28 


21576. 


.68 


146. 


.72 


2742. 


.55 


2208. 


.82 


178 


.18 


45 


.40 


33 


.26 
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Tabic 1: Masses of 32-dimensional even unimodular lattices (cont.) 



32256 
989231 

13440 
1223851 

46080 
6696269 
1182720 

2927 

2880 
31211779 
163296000 
1029287 
7464960 

66941 

184320 

3827 

8064 

34339 
384000 
783383 
13063680 
1021 

32256 
1 

252 

37 
2520 
4813 



dim root system R mW (R) ■ w(R) (as decimal) comments 

186.99 
73.60 
26.56 
5.66 
1.02 
0.19 
0.14 
0.36 
0.47 
0.09 
0.06 
0.03 
0.00 
0.01 
0.04 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 

0.00 unique 
0.00 
0.00 
0.00 
0.00 
0.00 

0.00 unique 
0.00 unique 
0.00 

0.00 unique 
0.00 shorter Leech 
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A A A'2 

A 1 A 2 Ai, 


9 


A a2 A 

A x AiA± 


9 


A2 A A 

AIA3A4, 


9 


A 9 A A 

A\A 2 A h 


9 


A\A\D± 


9 


A S 

A3 


9 


AfD 5 


9 


A 9 A T-\ 

A(A 3 D 4 


9 


A2A3A4 


9 


A A 

AiA 6 


9 


A 2 A 

AiA 5 


9 


A a2 


9 


A2A3D4 


9 


A\A 2 D 5 


9 


A1A3AT, 


9 


A!A 2 A 6 


9 


A1A4D4 


9 


AiD s 


9 


A 2 A 

AfA 7 


9 


A 1 A 3 D 5 


9 


A4A5 


9 


A 3 A 6 


9 


A 5 D 4 


9 


A 2 A 7 


9 


A{E 6 


9 


A 1 A 2 D 6 


9 


A4D5 


9 


AxA s 


9 


D 4 D 5 


9 


A 3 D 6 


9 


A\D, 


9 


A 3 E 6 


9 


A 9 


9 


A 2 E 7 


9 


D 9 



129024 

547 
126720 
1 

1260 
1 

8000 
43 
430080 
1 

15840 

13 
295680 
1 

466560 
1 

1512000 
1 

504000 
1 

1512000 
1 

645120 
1 

1774080 
1 

20401920 
1 

41287680 
1 

52254720 
1 

177408000 
1 

489646080 
1 

3592512000 

1 

991533312000 
1 

84610842624000 
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the order the automorphism group of the corresponding 26-dimensional lat- 
tice. Since mlj^CEe) • w(E 6 ) — |Aut(L )|~\ L must be the unique even 26- 
dimcnsional lattice of determinant 3 with no roots, and there must also be a 
unique 32-dimensional even unimodular lattice with root system Eq. (This also 
follows from the classification of 27-dimcnsional unimodular lattices with no 
roots in ||, as the one-to-one correspondence mentioned above also extends to 
27-dimensional unimodular lattices with no roots and with a parity vector of 
norm 3 j§.) 

Remark 8 In an earlier draft we pointed out the large mass of lattices with 
root systems A\ for small k, and remarked that roots seem to have a propensity 
for being orthogonal to one another. Peters has since sent us a preprint |}2| ex- 
plaining this: he observes that the sum of the masses m^ 2 (0) +Sfc=i m 32 C^i ) m 
| jl9| is 97.25% of the total mass of the genus of even unimodular 32-dimensional 
lattices, and shows that a lower bound of 97.11% of the total mass can be de- 
rived from the small size of the Fourier coefficient corresponding to A 2 in the 
Siegel Eisenstein series of degree 2 and weight 16. The idea is that an even 
unimodular lattice represents A 2 if and only if it has non-orthogonal roots u 
and v (with u ^ ±v). 

4 Masses of unimodular lattices of dimension 
n < 30 with any given root system 

In ||, Chapter 16], Conway and Sloane describe a correspondence between uni- 
modular lattices of dimension n < 23 and orbits of norm 4 vectors in even 
unimodular 24-dimensional lattices. They use this correspondence to produce, 
from the list of Niemeier lattices, a list of the root systems and automorphism 
group orders of all unimodular lattices of dimension n < 23. Unimodular lat- 
tices of dimension n < 31 likewise correspond to orbits of norm 4 vectors in 
even unimodular 32-dimensional lattices, and we can use this correspondence to 
compute m n (R) for all root systems R and all n < 30. 

Let v = 2e be a vector of norm 4 in an even unimodular lattice A of dimension 
n = 32. Then L31 = {x E A 1 - : x + ne E A for some n E Z} is an odd 31- 
dimensional unimodular lattice, and |Aut(L3i)| = 2|Aut(A)|/c(w), where c(v) is 
the number of images of v under Aut(A). If the lattice L31 has exactly 2k vectors 
of norm 1, we can write L31 = L%i^ r © Z fc , where L^i^k has minimal norm 2 
and |Aut(L 3 i_fc)| = |Aut(L 3 i)|/(2 fc fc!). (We shall call L^i-k the reduced lattice 
corresponding to A and v.) This construction gives a one-to-one correspondence 
between orbits of norm 4 vectors in 32-dimensional even unimodular lattices and 
unimodular lattices with no vectors of norm 1 in dimensions less than 32. (See 
[gor§). 

We shall mainly be concerned with norm 4 vectors v which are the sum of 
two orthogonal roots r and s. For such vectors, knowing just the root system of 
A allows us to compute the dimension and the root system of the reduced lattice 
Lai—k- Some useful information, distilled from J9|, Table 16.8], is provided in 
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the table below. 

R #r R shape of v #w R dim(L 3 i_fc) 



A n 




An-2 


any 


err) 


A n — 4 


29 


D 4 


24 


A\ 


any 


24 





28 


D n (n > 5) 


4© 


Da-2 


(±i 4 ,o n - 4 ) 

(±2,0"- 1 ) 


16Q) 
2n 


D n -i 



28 
32 


Eq 


72 


A 5 


any 


270 





27 


E 7 


126 


D e 


any 


756 




26 




240 


E 7 


any 


2160 





24 



The interpretation of this table is as follows: 

Case 1: r and s are from different components R and S of the root system 
of A. Then the reduced lattice L^\-k has dimension 30, and has the same root 
system as A, except for the components R and S from which r and s are taken, 
which are transformed to R and S as given in the third column. The second 
column gives the number of roots #r in each component R, from which we can 
compute the number of vectors v — r + s with r e R and s e S. 

Case 2: r and s are from the same component R of the root system of 
A. Then the dimension of the reduced lattice is given in last column, and the 
root system of L^-k is the same as that of A except the component R will be 
replaced by R. Note that for R = D n with n > 5, R and dim(L 32 _fc) depend 
on the shape of v. The column headed #w gives the number of norm 4 vectors 
of that shape in the component R. 

Remark 9 Norm 4 vectors that are not the sum of two orthogonal roots corre- 
spond to reduced lattices of dimension 31. For this reason we have not computed 
the mass of unimodular 31-dimensional lattices with root system R. (That mass 
could be computed from the mass of even unimodular 40-dimensional lattices 
with root system R © Dg.) 

Computing the number c(v) (and hence the order of Aut(Zv3i_fc) in terms of 
the order of Aut(A)) sometimes requires additional information about A, such 
as whether there are automorphisms of A which permute the components of its 
root system with multiplicity greater than one. But for our purposes this is not 
an impediment, as we demonstrate in the following example. 

Example 10 Let A be an even unimodular 32-dimensional lattice with root sys- 
tem AfD$ and let v = s+t where s is a root from one of the four components A\ 
and Ms a root from the component D 5 . Then the corresponding reduced lattice 
L- 3 i-k is 30-dimensional, has root system AfA 3 (since D 3 — As), and has auto- 
morphism group order |Aut(i 30 )| = (2|Aut(A)|/c(w))/(2 1 • l!) = |Aut(A)|/c(u). 
Since each component A\ has 2 roots and the component D 5 has 40 roots, there 
are 4 • 2 • 40 = 320 such vectors v, which form anywhere from one to four orbits 
under Aut(A), depending on whether there are any automorphisms permuting 
the components A\. Because of this ambiguity, we do not know exactly what 
c(v) is. But suppose the 320 vectors break into m orbits, Vi,...,V m , with 
representatives v\, . . . , v m . Then 

m m 

]t C (^) = En = 320 

i=l i=l 
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so the total mass of the lattices which correspond to the lattice A and any of 
these vectors v is 



Since from Table [j] the mass of all 32-dimensional lattices Aj with root sys- 
tem AfD 5 is J2j lAutCAjOr 1 = 1029287/7464960, these lattices contribute 
£V 320|Aut(A J )|- 1 = 320 • 1029287/7464960 towards the total mass of 30- 
dimcnsional unimodular lattices with root system A3. By similarly account- 
ing for contributions from orbits of norm 4 vectors in even unimodular 32- 
dimensional lattices with other root systems, we can compute the exact mass of 
the 30-dimensional unimodular lattices with root system 

We have, in this manner, computed the mass of n-dimensional unimodular 
lattices having any given root system for all n < 30. We have not provided a 
table of these masses (as they can easily be derived from the table of masses 
of 32-dimensional even unimodular lattices with any given root system in fl9[| ), 
but we will use these masses in Section || to find lower bounds on the class 
numbers of unimodular lattices in dimensions up to 30. 

Example 11 Even unimodular 24-dimensional lattices with root system R cor- 
respond to even unimodular 32-dimensional lattices with root system R © Eg, 
and m 2 4 (R) = m 3 \ (R®E$) -w(E$) -c where c is the multiplicity of E$ in R®Eg. 
We see from the table in [jl9| that there are 24 root systems of the form R © Eg 
for which m 32 {R®E 8 ) > 0. Since even unimodular 24-dimensional lattices hap- 
pen to be uniquely determined by their root systems, the corresponding values 
of m.24 (-R) are precisely of the form |Aut(A)| _1 , where A is the Niemeier lattice 
having root system R. 

5 Mass formula for unimodular lattices with no 
roots 

Recall that m„(0) denotes the mass of the n-dimensional unimodular lattices 
with no roots. By considering all the orbits of norm 4 vectors v = r + s in even 
unimodular 32-dimcnsional lattices that correspond to reduced lattices with no 
roots, we have 

m n (0) = m\ l 2 {D 32 - n )w{D 32 - n ) for n < 26, n ^ 24 



m 24 (0) 


= m L 3 l 2 {D 8 )w{D 8 )^ 


-m 3 l 2 (E 8 )w{E 8 ) 


m 27 (0) 


= m l i 2 {D 5 )w{D 5 )^ 


-mg(E e )w(E e ) 


m 28 (0) 


= m 3 I 2 (D 4 )3w(Z? 4 ) 


+ m\ 1 2 (D 5 )w{D 5 ) 


m 29 (0) 


= mf 2 (A 3 )w(A 3 ) + 


77132 (At) u>(Ai) 


m 3O (0) 


= rr$ 2 (Al)w(Al) + 


m\\ {A X A 2 ) w{A 1 A 2 ) + m\\ (A\) w{A 2 2 ) 



These masses m„(0) may be computed by looking up the values of m\\{R) 
in Table [l] and p9[ ; the results are listed in Table ||, split into 771^(0) and m^ 1 (0) 
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Table 2: Masses of n-dimensional unimodular lattices without roots. Odd lat- 
tices appear above, even lattices below. 



(as decimal) number of lattices 



0-22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 



84610842624000 

1 

1002795171840 



1 

18720000 
206867 



1585059840 
17924389897 
26202009600 
49612728929 
11136000 
7180069576834562839 
175111372800 

? 
? 
? 








-14 


1 


-13 


1 







-8 


1 


-4 


3 


-1 


38 



4.46 x 10 3 
4.10 x 10 7 

> 4 x 10 11 

> 5 x 10 15 

> 4 x 10 20 



(shorter Leech lattice) 
(odd Leech lattice) 

(classified in |I|) 
(classified in j2|) 
(classified in |||) 
more than 8900 
more than 82000000 
more than 8 x 10 11 
more than 1 x 10 16 
more than 8 x 10 20 [i 



(as decimal) number of lattices 



16 
24 
32 



8315553613086720000 
1310037331282023326658917 
238863431761920000 



1 



1.20 x 10~ 19 
5.48 x 10 6 



1 (empty lattice) 





1 (Leech lattice) 
more than 10000000 
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for n divisible by 8. Table ^ also includes lower bounds on m l n {$) for n = 31 
and 32, which we explain below. 

Suppose A is an even unimodular 32-dimcnsional lattice with no roots. It 
follows from a theta function argument as in ||, Theorem 7.17] that A has 146880 
vectors v of norm 4. For each orbit of these vectors, the corresponding reduced 
lattice is 31-dimensional and has no roots. Thus these lattices A contribute 

146880 TT ... 22270634631794396553201589 „ „„ 

m" (0) = w 4.03 x 10 11 

2 32 w 55292461056000 

towards the mass of 31-dimensional lattices with no roots. There are almost 
certainly additional contributions from norm 4 glue vectors v in even unimodular 
32-dimensional lattices with root systems A\ for some k, but we will not attempt 
to account for these. We can restate this in terms of parity vectors: 



Proposition 12 The mass of SI -dimensional unimodular lattices with no roots 

4W- 



and with no parity vectors of norm 7 is (146880/2) 



Proof. Let L be a 31-dimensional unimodular lattice with no roots, and A the 
corresponding 32-dimensional even unimodular lattice. As in ||, p. 414], A is 
equal to {x + y \ x G L\ y G \ + 2Z, for i = 0, 1, 2, 3} C (L )' ±Z, where L° 
is the sub-lattice of L consisting of vectors of even norm; L°, L 1 , L? and L 3 are 
the cosets of L° in its dual (L )'; and L° U L 2 = L. The only ways the vector 
x + y in A can have norm 2 is if y = and (x, x) — 2 with x G L , if y = 1 and 
(x, x) = 1 with x G L 2 , or if y = ±1/2 and (x, x) = 7/4 with x G L 1 or x G L 3 . 
In the first two cases x would be in L which would contradict L having no roots. 
Thus A has a root if and only if there is a vector x of norm 7/4 in L 1 U L 3 . But 
in this case x is in {LP)' \ L — the so-called shadow of L — and 2x has norm 7. 
Since the parity vectors of L are precisely twice the shadow vectors |], Preface 
to 3rd edition, p. xxxiv], L has a parity vector of norm 7 if and only if A has 
a root. (The 31-dimensional unimodular lattices with no roots and no parity 



vectors of norm 7 must have parity vectors of norm 15 by 1 1 1 .) 



We can construct 32-dimensional odd unimodular lattices with no roots from 
32-dimensional even unimodular lattices with no roots as follows: Let A be a 
32-dimensional even unimodular lattice with no roots, and let b € A/2A be a 
nonzero element with norm divisible by 4. Then there is a unique 32-dimensional 
odd lattice L containing A b := {v G A | (v,b) G 2Z} (see [§J Chapter 0.2]). (A 
and L are neighbors, meaning their intersection has index two in each of them.) 
If b is not represented by a vector of norm or 4 then L has no roots. Of the 
2 32 elements in A/2A, 2 31 + 2 15 have norm congruent to mod 4, by Milgram's 
formula |2^, Appendix 4] applied to the lattice \J~2 A. (Milgram's formula says 
that if M is an even lattice and M' its dual, then 

^2 cxp(27ri(u,u)/2) = (det Af) 1/2 exp(27ricr/8) 

ueM'/M 

where a is the signature of M, which is equal to its dimension when M is 
positive definite.) Note that (6,6) mod 4 depends only on 6 mod 2A. Vectors 
v and w of norm 4 in A are equivalent mod 2A if and only if v = ±uj, so 
146880/2 elements of A/2A are represented by vectors of norm 4. Thus there 
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are 2 31 + 2 15 - 146880/2 - 1 = 2147442975 elements of A/2A represented by 
vectors of norm divisible by 4 but not represented by vectors of norm or 4. 
Since each odd 32-dimcnsional unimodular lattice has two even neighbors, this 
implies that m\ 2 {%) > (2147442975/2) • m| I 2 (0) = 5.89 x 10 15 . (This is almost 
certainly an underestimate, since an even neighbor of an odd lattice with no 
roots can have roots.) 

As an immediate consequence of the values we computed for m„(0) in Ta- 
ble ||, we have: 

Corollary 13 ([^, There exist odd unimodular n- dimensional lattices 

without roots for n = 23, 24, and 26 — 32, but not for n = 1 — 22 or n = 25. 

Unimodular lattices with no roots are known to exist in all dimensions n > 
26 (see B, [29 ). In dimensions n < 28 they have already been completely 



enumerated, and so the masses m„(0) may also be computed by summing the 
reciprocals of the automorphism group orders of these lattices. (Happily, this 
agrees with our mass formula in each case.) For the even lattices, there is the 
empty lattice in dimension 0, and the Leech lattice A24 in dimension 24. The 
Leech lattice A24 was discovered by Leech in 1965 |^6| and was shown to be the 
unique even unimodular 24-dimensional lattice without roots by Niemeier |50| 
and by Conway Chapter 12] around 1969. For the odd lattices, there is the 
shorter Leech lattice in dimension 23, the odd Leech lattice O24 in dimension 24 
[ pl| , the lattice S26 in dimension 26, 3 lattices in dimension 27, and 38 lattices 
in dimension 28. The lattice 5*26 was constructed by Conway in the 1970's and 
was shown to be the unique unimodular 26-dimensional lattice without roots by 
Borcherds in 1984 Q|. Borcherds also found one of the 27-dimensional lattices. 
The full enumerations in dimensions 27 and 28 are due to Bacher and Venkov 



Remark 14 In the cases of the Leech lattice and S26, with our mass formula 
the uniqueness follows immediately from the constructions, simply by verifying 
that |Aut(A 24 )r 1 = m| T 4 (0) and lAut^e)! -1 = m \&{®)- 

Examples of unimodular lattices with no roots have been constructed for di- 
mensions 29 to 32 (including the 15 exceptional even unimodular 32-dimensional 
lattices classified in [p5| ), and a nonconstructive analytic argument shows that 
they exist in all dimensions n > 33 (see the Conway-Thompson Theorem p8[ 
p. 46] for n > 37, and Q for 33 < n < 36.) In fact for n > 33 this argument 
gives a lower bound for the mass of lattices without roots that is close to the 
total mass of the genus, so there are a great many lattices without roots (see 
Remark |l~6|). The idea is that the coefficients a\ and 02 of the average theta 
series 

- TT"77TTT = 1 + a ^ + a ^ + a ^ + '~ 
m tti ' Aut ( A )l 

give the average number of vectors of norm 1 and norm 2, taken over all the 
lattices in the genus f2. If a\ + a-i is less than 2, then there must be some A e SI 
with no vectors of norm 1 or 2 — this is the case when is the genus of n- 
dimensional unimodular lattices for n > 33. For n — 33, 01+02 is approximately 
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1.42, and this implies that m\ 3 ($) > 4.04 x 10 20 , so there are more than 8 x 10 20 
33-dimensional unimodular lattices without roots For n < 32, the average 
number of roots is greater than 2, so this argument does not apply. Notice that 
the coefficient Oi is the same as the average number of representations a(N) we 
defined in Section ^, in the special case where N is the lxl matrix (i). 

Remark 15 In dimensions n < 33, odd unimodular lattices with no roots have 
minimal norm 3, except in dimension 32, in which they can have minimal norm 
3 or 4. In dimensions 24 and 32, even unimodular lattices with no roots have 
minimal norm 4. See Hi. 



Remark 16 Let m' n be the mass of odd n-dimensional unimodular lattices 
with only trivial automorphisms. There are no such lattices for n < 28 but 
Bacher has found one for n = 29 0. Bannai || showed that m' n /m n — > 1 as 
n — > oo. For n > 1 any lattice with roots has nontrivial automorphisms, so 
m n/ m " — m n(0)/^ri- Below we list to„(0)/to„ for 26 < n < 30, and lower 
bounds on m n (0)/m„ for 31 < n < 33. 

n m„(0)/m„ ref. n m„(0)/m„ ref. 



26 0.000116 \4] 30 0.0908 

27 0.000856 2] 31 > 0.135 

28 0.00658 § 32 > 0.136 

29 0.0300 33 > 0.287 § 

6 Lower bounds on class numbers 

Let Q be the set of inequivalent lattices in a genus of dimension n > 0, m the 
mass of that genus, and m{R) the mass of those lattices having root system R. 
Each lattice A € has at least two automorphisms, 1 : x i— * x and —1 : x i— ► — x; 
from this we get the well-known lower bound |0| > [2m] . For each root r of A, 
the reflection 

x h- > x — 2- -r 

{r, r) 

is also in Aut(A). Define w'(R) to be the order of the subgroup of Aut(A) 
generated by reflections and by — 1. Then w'(R) = w(R) if the map —1 is 
already in the Weyl group of R, and w'(R) = 2w(R) otherwise. (The Weyl 
group of R contains —1 if and only if rank(i?) = dim (A) and each component of 
R is Ai,Et, Eg or Dk for even k.) Then there are at least \m(R) w'(R)~\ lattices 
with root system R, so we get an improved lower bound, 

\n\>J2\m(R)w'(R)-}. 

R 

We can do slightly better still, as follows: Write m(R) w'(R) = q + a/b with 
5, a, b G Z, a < 6, and gcd(a,6) = 1, and define a modified ceiling function by 
(q + a/b) = q if a = 0, q + 1 if a = 1, and q + 2 if a > 1. It can easily be shown 
that there are at least (m(R) w'(R)) lattices with root system R, so that 

M>J2{m(R)w'(R)). (1) 
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Evaluating this sum with the value of m l ^ 2 (R) we computed for each R gives 
|fi| > 1162109024. 

Corollary 17 There are at least 1000000000 even unimodular 32 -dimensional 
lattices. 

For comparison, [277132] = 80618466, so our lower bound is 14.4 times larger 
than the lower bound obtained by doubling the Minkowski-Siegel mass constant. 

We also computed m n (R) for each n < 30 and each root system R, and used 
equation [l] to find lower bounds on the number of odd unimodular lattices. (In 
the cases where a single m„(i?) is expressed as a sum of masses corresponding 
to different root systems of 32-dimcnsional even unimodular lattices, such as 
TO2s(0) = W32 (-D4) 3u>(Z?4) + m 3 2 (D5) w(Ds), we bounded the number for each 
summand individualy, which gives a better overall bound.) 

Unimodular lattices in dimensions n < 25 have been completely enumerated: 
even unimodular lattices of dimension 8 by Mordell, of dimension 16 by Witt, 
and of dimension 24 by Niemeier ||(J (see also Venkov ||3]); odd unimodular 
lattices of dimension n < 16 by Kneser pi] , of dimension n < 23 by Conway 
and Sloane @, Chapter 16], and of dimension 24 and 25 by Borcherds Q. 

Table | gives our lower bound n on the number of unimodular lattices in 
dimension n < 30, and our computation of the number of distinct root systems 
r n that occur in these lattices (including root systems with components Z, which 
do not occur for even lattices). The table includes for comparison the actual 
number a n of unimodular lattices in dimension n < 25, taken from §, Table 2.2], 
and the Minkowski-Siegel mass constants m n , taken from || Tables 16.3 and 
16.5]. (The counts of odd unimodular lattices include those with vectors of norm 
1 to facilitate comparison with m n ; since any integral lattice with a vector of 
norm 1 is of the form Z fc © A where A has minimal (nonzero) norm 2, the counts 
of lattices with no vectors of norm 1 can be recovered from /3„; similarly, the 
number of distinct root systems with no components Z can be recovered from r n , 
with the caution that 8 of the 24 root systems of even unimodular 24-dimensional 
lattices also occur as root systems of odd unimodular 24-dimcnsional lattices.) 
Our lower bounds agree exactly with the actual numbers for n < 24. This 
is to be expected for n < 23 since in those dimensions an odd or even n- 
dimcnsional unimodular lattice is uniquely determined by its root system i?, 
and {m(R)w'(R)} is exactly equal to the number of lattices with root system 
R when there are or 1 such lattices. This is also true for 24-dimensional even 
unimodular lattices. Our lower bound in dimension 25 is within two percent of 
the actual number, and our lower bounds in dimensions 26 to 30 (for which the 
actual numbers are not known) are the best we are aware of. 

7 Computing the numbers a(R) 

For a half- integral n x n matrix B, define 

c„. fc (S) = (-l)" fc/2 2"( fc -^ 1 )/ 2 )(det73)( 2fc — 1 )/ 2 &(S,A : ) TT 

- L - L 1 (1/2) 



17 



Table 3: Comparison of the actual number a n of unimodular lattices of di- 
mension n, the lower bound (3 n computed with equation |], the number r n of 
distinct root systems occuring (including those with components Z), and the 
Minkowski-Siegel mass constant m n . Odd lattices are listed above, even lattices 
below. The last column is included only when 2m n provides a nontrivial lower 
bound for a n . 
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where 

b(B,k)= exp(2mtT{BR))fi(R)- k 
fles n (Q)/s n (z) 

is the Siegel series (with equal to the product of denominators of ele- 

mentary divisors of R). Put e„.k = 1/2 if n = k — 1 or n = k > 1, and 1 
otherwise. 

Theorem 18 (Siegel) Let N be a lattice with dim(TV) = n < 8fe, and ie< 



r(A, JV) 



< aVo > Aut ( A )l 

imt/i £fte sum taken over even unimodular lattices of dimension 8k. Then a(N) = 
£n,skC n ,4k(B) where B is 1/2 times the Gram matrix of N . 

Proof. See Theorem 6.8.1], and note that the product of local densities 
Yl„ QipjMp, Np) is equal to b(M,k) when N is an even unimodular lattice of 
dimension 2k |2l[| . 

Remark 19 The Siegel Eisenstein series of degree n and weight k is defined 
to be 

E n , k (Z)= \CZ + D\- k 

{C,D} 

where {C, D} runs over all representatives of the equivalence classes of coprime 
pairs of n x n matrices. For k > n, c n ,k{C) is the coefficient in the Fourier 
expansion 

E n , k {Z) = ^c n , fc (C , )exp(27ritr(CZ)) 
c 

where C runs over all positive semi-definite half-integral n x n matrices. 

The first explicit formula for b(B, k) for arbitrary n is due to Katsurada ||i6f , 
and was published in 1999. Prior to this, an explicit formula for the coefficients 
of E n ,k(Z) was known only for n < 3: The case n = 1 is well-known (see for 
example jB3[ Chapter VII]), Maafi 27 gives an explicit formula for n = 2 (see 
also Kaufhold [l7|), and Katsurada 



(extending partial results by Kitaoka 



[ lq| gives an explicit formula for n = 3 

Let _B be a nondegenerate symmetric half- integral n x n matrix over Z p . It 
follows from |E3| that 



« B ,.)=(<(.) l nc( 2 .- 2i) ) n ^->{f ( *" r/2;M) 

\ i=l / P\D(B) K - 

for certain polynomials F p (B;X). Here C is the Riemann zeta function, and 
L(s, xb) is a Dirichlet L-series whose values may be computed using the method 
in §. 

The explicit formula for F P (B;X) in |l(| Theorem 4.3] is not itself well- 
suited for calculation since the outer index of summation takes 2™ values, but 
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we can use Katsurada's recursion relations Jl6|, Theorems 4.1 and 4.2] as part 
of a practical algorithm for computing F p (B;X). We will state these recursion 
relations below (without proof); first we will need to introduce some of the 
notation from pq |. 

For a = p r c with r € Z and c € Z*, define Xp( a ) — (§) f° r r even and for 
r odd (where (-) is the Legendre symbol mod p), and define 

{+1 if r = mod 2, c = 1 mod 8 
— 1 if r = mod 2, c = 5 mod 8 
otherwise. 

Define ord p (a) to be the exact power of p dividing a, and define i p (B) to be 
the least integer t for which p t B^ 1 is half-integral. Let ( , ) p denote the Hilbert 
Symbol over Q p (see |53|), and let h p denote the Hasse invariant (see p2|). For 
odd n, define 

Vp (B) = h p (B)(dctB,(-l) {n - 1)/2 detB) p (-l,-l)f^/ s . 
For even n, define 

UB)=X P ((-l) n/2 detB) 

and C P (B) = 1 + C P {B) - £ P (B) 2 . By convention, £ p (B) = £' p (B) = 1 for B the 
empty matrix. Define D(B) = 2 2 L"/2J detB, d p (B) = oid p (D(B)), and 



S P (B) 



2[(d p (B) + 1 - 5 2p )/2\ if n is even 
d p (B) if n is odd 



where 52 P is the Kronecker delta. 

We shall suppress most of the subscripts p in what follows. Let p be any 
prime and suppose B and B 2 are nondegenerate half integral matrices of rank 
n and n — 1 respectively over Z p . Put <5 = 5(B) and <5 = 5(B 2 ). If ft* is even, put 
£ - £(5), £' - ?(B), and ry = t?^); if n is odd, put £ = £(B 2 ), |> = £'(B 2 ), 
and 77 = T)(B). (By convention 77 = 77' = 1 if _B 2 = 0, and 5 = if n = 1.) Then 
define rational functions C(B, B 2 ; X)^ and C(B, B 2 ; X) (0 *> in X by 

c(5,s 2 ;x)« = { ^ isev ;; 

1 11 71 is odd 



l_p(«+l)/2£ X 



and 



ifniseven 

( { { _^ (v{n _^^ p{ ^ )/2 . fnisodd 



Theorem 20 (Katsurada ||l6| ) Let £?i = (61) and B 2 he nondegenerate half 
integral matrices of degree 1 and n — 1, respectively, over Z p , and put B = 
£?i_L£? 2 . Assume that ord(b\) > i(B 2 ) — 1 + 25 2p - TTien we /lave 

F P (B; X) = C(B, B 2 ;X)^F(B 2 ;pX) + C(B, B 2 ; X)<® F p (B 2 ; X). 
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Let B 2 be a nondegenerate half-integral matrix of degree n — 2 over Z 2 , let 
H = (° 1/2 1/2 ) and let Y = ({^l)- Let B t = 2 m K with K = H or F, or ^ = 

2 m u 1 _L2 m u 2 with u x ,it 2 S Z£, and put 5 = Bi_LB 2 , 5 = 5(B), 5 = 8(2 m ±B 2 ), 
5 = 5{B 2 ), and 

(25 - S - 5 + 2)/2 if n is even, Bi = 2 m u 1 _L2 m w 2 and d(B) is odd 

or if n is even, B x = 2 m K, and £(B 2 ) = 
2 if n is odd, B x = 2 rn K, and d(2 m ±B 2 ) is even 

otherwise. 

If n is even, put £ = £(B), f - £'(£), £ = £(#2), I' = £'(£? 2 ), and 

' r]{2 m u 2 ±B 2 ) if £1 = 2 m wi_L2 m w 2 , 

and d(B 2 ) is even 
?? = <| (-l)«"- 1 ) 2 - 1 )/ 8 /i(S 2 )(2 m , (-l)^- 1 )/ 2 det B 2 ) 2 if B 1 = 2 m K, 

and £(B 2 ) ^ 

1 otherwise. 

If n is odd, put 7/ = rj(B), i) = i]{B 2 ), £' = 1, and f = 1 if B x = 2 m K and 
d(2 rn -LB 2 ) is even, and otherwise. Define four rational functions in X by 



C(B,B 2 ;X)^ = 
C(B,B 2 -X)^ = 



1-2-WP ^ U 1S even 

1 _ 2( , + 1 i)/2; X if "is odd 



l- 2 n+lX 2 

l- 2 ("+l)/2^X 



if n is even 
if n is odd 



1 



. if n is even 

C(B,B 2 ;X) {21 > = < ll 2 "«-D/2j X . 

I 1-2"^ — 71 1S 

(_ 1 * ) £{:lf l ( 2 (n-2)/2 x y-5 + 2-i 2 -<, 2 (25-S + 2-2<T)/2 



if n is even 



C(B,B 2 -X)^ = { _ i-a»/»£x . 

1 (-i)^^(i- 2 '" +1 ^^)g : -fx)-^ 2 - 2 '^^) if n is odd . 

Theorem 21 (Katsurada |l§) Let B x = 2 m u 1 ±2 m u 2 with m,u 2 E Z* 2 or 
B = 2 m K with K = H or Y . Let B 2 be a half-integral matrix of degree n — 2 
over Z 2 which is also in Gi„_ 2 (Q 2 ), and put B = B\JlB 2 . Assume that m > 
i(B 2 ) + 1. Then we have 

F 2 (B;X) = C(B,B 2 ;X)^C(B,B 2 ;2X)^F 2 (B 2 ;4:X) 

+ C(B, B 2 ; X)^C{B, B 2 ; 2X)^F 2 (B 2 ; 2X) 

+ C{B, B 2 ; X)WC{B, B 2 -X)^F 2 {B 2 - 2X) 

+ C(B, B 2] X)( W) C{B, B 2 -X)^F 2 {B 2 ;X). 

For p 7^ 2, any nondegenerate symmetric half-integral nx n matrix B can be 
diagonalized over the p-adic integers Z p , B = p ei u\±- . . . ±p e "u n with e\ > ■ ■ ■ > 
e n > and with Ui S {l,e} for all i, where e is any quadratic nonresidue (see 
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[p5|). Then F p (B; X) can be computed by repeated applications of Theorem p0| 
above. 

Note that in the course of recursively computing F p (B;X) for X = p~ k 
in this manner, one occasionally encounters zeros of the denominators of the 
functions and . (This is only a problem when X is a negative power 
of p.) Rather than attempting to simplify the expressions symbolically we 
instead computed F p (B; 1), F p (B;p), . . . , F p (B;p n ), and then used Lagrangian 
interpolation to compute F p (B;p~ k ). The reason for choosing X = l,p, . . . ,p n 
rather than n + 1 other numbers is that for these numbers the recursive sub- 
problems overlap: we need only evaluate F p (p ei Ui-L . . . -Lp en u n ; X) for X = 
l,p, . . . , p 2 ™ -4 , starting with i — n, and working down to i = 1 which gives 
F p (B;l),...,F p (B;p n ). 

For p = 2, any nondegenerate symmetric half-integral matrix B is equivalent 

over Z 2 to a matrix of the form 2 ei (Ui±Vi)J L2 e ™ (U m ±V m ) with e\ > ■ ■ ■ > 

e m > 0, Ui = or u\ or ui±u 2 for ui, u% 6 {±1, ±3}, and Vi = or H± ■ ■ ■ ±H 
or H± ■ ■ ■ _LiJ_l_y (see p5|]). Then F 2 (B;X) can be computed by repeated 
applications of Theorems 20 and |2l| above. (We again used an interpolation 
scheme, similar to the one described for p ^ 2.) 



8 Computing the numbers r(Ri, Rj) 

The method described in this section is essentially the one used for some of 
the computations in || (although the algorithm itself is not described in that 
paper). In this section we write r(R,R') as emb(R',R), since a representation 
of R 1 by R is the same as a linear map from R' into R which preserves inner 
products. For any irreducible root systems S and T it is routine to compute the 
number of embeddings emb(5, T) of S into T, and to determine the root system 
of the orthogonal complement of S in T for each of these embeddings. There 
are at most two orbits of embeddings of S into T, so we write emb(5, T) = 
embi(S l , T) + emb2(S,T), where there are embi(5, T) embeddings of S into T 
for which the orthogonal complement is comp 1 (S r , T) and there are enu^S 1 , T) 
embeddings of S into T for which the orthogonal complement is comp 2 (S r , T). 

Example 22 We demonstrate how to compute the two orbits of embeddings of 
A3 into D n . Recall from H that A n = {(xi, . . . , x n+ i) e Z™ +1 : x\ + - ■ -+x n+ i = 

0} and D n — {{x\, . . . , x n ) £ Z n ; x% H + x n e 2Z}. A 3 is generated by 

three roots v x = (1,-1,0,0), v 2 = (0,1,-1,0), and v 3 = (0,0,1,-1) with 
(vi,vi) = (v 2 ,v 2 ) = (v 3 ,v 3 ) = 2, (vi,v 2 ) = {v2,v 3 ) = -1 and («i,u 3 ) = 0. 
We can map v\ to any of the 4C(n, 2) roots r% in D n (all permutations of 
(±1, ±1, 0, . . . , 0)), where C(n, k) denotes the binomial coefficient. 

We can then map v 2 to any of the 4(n — 2) roots r 2 of D n which have inner 
product —1 with r\. Let us say T\ is supported in coordinates i and j and r 2 is 
supported in coordinates j and k. (Clearly k ^ i.) Then there are two cases: 

(a) We can map i>3 to the root r 3 supported in coordinates i and j which has 
inner product with n and inner product —1 with r 2 . In this case the roots of 
D n orthogonal to r±, r 2 , and r 3 form the system D n _ 3 ; 

(b) We can map v 3 to 2(n — 3) roots r 3 supported in coordinates k and 
/ (I which have inner product with n and inner product —1 with r 2 . 
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In this case the roots of D n orthogonal to n, r%, and r% form the system Z) n _4. 

Hence there are 4C(n, 2) ■ 4(n — 2) • 1 = C(n, 3) ■ |Aut(^3)| ways to embed 
A3 into D n with complement -D n -3, and there are 4C(n, 2) • 4(n— 2) • 2(n— 3) = 
2 3 C(n,4) • |Aut(A 3 )| ways to embed A 3 into D„ with complement D„_4. 

The computations for other irreducible root systems are similar. All nonzero 
values of embi, emb2, comp^ and comp 2 are given in Table |]. 

If R and R 1 are root lattices then we can write R = S% © • • ■ © Sk and 
R' = Tx © • • • © T m where each S 1 , and is an irreducible root lattice, and we 
can compute r{R' , R) recursively via the formula 

r(R',R) = emb(R,R') 

= emb(05 i ,0T i ) 

i<k £<m 

rn 

= ^emb 1 (5 t ,T J )emb(0 S^comp^St,^)® I}) 

j — l i<k — l i<mj,^j 

rn 

+ ^cmb 2 (S fe ,r j )emb( comp 2 (S fe , 3}) © Tj). 

j — l i<k — 1 i<m,i^j 

A direct implementation of this algorithm does a lot of redundant compu- 
tation on certain inputs, some of which we can circumvent with dynamic pro- 
gramming or memoization (sec |l0| , Chapter 16]). Our implementation also does 
several things to reduce the amount of computation when there are direct sum- 
mands in R or R' with multiplicity greater than one. But computing r(R' , R) 
is an NP-Hard problem, since it is an NP-Complete problem to determine if 
r{R' , R) > 0, as we show below. 

Proposition 23 The problem of determining whether a root system R embeds 
into a root system R' is NP-Complete. 

Proof. We shall reduce 3-Partition (see @) to this problem. Let S — 
{si, . . . , ssk} be an instance of 3-Partition, with the s$ positive integers sum- 
ming to kt. Then S can be partitioned into k sets each consisting of 3 elements 
with sum t if and only if the root system D± Sl © • • • © D^ S3k embeds into the 
root system © i=1 D^. 

Remark 24 Since 3-Partition is strongly NP-Complete, the problem of whe- 
ther one root system embeds into another remains NP-Complete if the Gram 
matrices of root lattices are used as input (rather that the list of components 
Ai, Di, and Ei as above). 

Since we need to compute r{R, R') for all pairs of root systems rather than 
just one pair, the amortized computational cost would be reduced considerably 
by using dynamic programming. But the dynamic programming table becomes 
unmanageably large in dimension 32, so we instead use a hash table, which is 
purged periodically, for memoization. 
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Table 4: Embeddings of irreducible root systems into one another other. (In 
the fourth and last columns, A , D , D\, D 2 and D 3 should be interpreted as 



0,0,0,^1^1 and A 3 , respectively.) 



b 


1 


cmbi(S.T) 
|Aut(S)| 


comp^o, 1 ) 


cmb 2 (S,T) 
|Aut(S)| 


comp 2 (»b,i ) 





T 


1 


T 


— 


— 


M 


Aj (j > i) 


C(j + 1, i + 1) 


Aj-i-i 





— 


Ai 


Dj 


C(j, 2) • 2 


AiDj-2 


— 


— 


A 3 


Dj 


C(j,4) -2 3 


Dj-4 


C(j, 3) 


Dj-3 


Ai (i 7^ 1,3) 


Dj (j > i) 


C(j, i + l)-2' 


Dj-i-i 


— 


— 


A 


Dj (j > i) 


C(j,i) 

_ 9 „ 9 

2 2 3 


Dj-i 


— 


— 


Ai 


Ee 


A r , 


— 


— 


A 2 


Ee 


2 3 3 • 5 


A 2 A 2 


— 


— 


A 3 


Ee 


2 ■ 3 3 5 




— 


— 


Ai 


Ee 


2 3 3 


A! 


— 


— 


A 5 


Ee 


2 3 


A! 


— 


— 


D A 


Ee 


3 2 5 





— 


— 


D 5 


Ee 


3 3 





— 


— 


Ee 


Ee 


1 





— 


— 


Ai 


E 7 


3 2 7 


D 6 


— 


— 


A 2 


E 7 


2 4 3 • 7 


A,, 


— 


— 


A 3 


E 7 


2 2 3 2 5 • 7 


A 3 Ai 


— 


— 


A A 


E 7 


2 5 3 2 7 


A 2 


— 


— 


As 


E 7 


2 4 3 • 7 


A 2 


2 4 3 2 7 


Ai 


A 6 


E 7 


2 5 3 2 





— 


— 


A 7 


E 7 


^.9 „9 

2 2 3 2 





— 


— 


D 4 


E 7 


3 2 5 • 7 


AiAiAi 


— 


— 


D 5 


E 7 


2 • 3 3 7 


Ai 


— 


— 


D 6 


E 7 


„ 9 _ 

3 2 7 


A! 


— 


— 


Ee 


E 7 


~9 _ 

2 2 7 





— 


— 


E 7 


E 7 


1 





— 


— 


Ai 


Eg 


2 3 3 • 5 


Ei 


— 


— 


A 2 


E& 


2 5 5 • 7 


Eq 






A z 


Eg 


2 3 3 3 5 ■ 7 


D 5 






A 4 


Eg 


2 7 3 3 7 


A 4 






As 


Eg 


2 7 3 2 5 • 7 


A 2 Ai 






A 6 


Eg 


2 8 3 3 5 


A ± 






A 7 


Eg 


2 5 3 3 5 


Ai 


2 6 3 3 5 





As 


Eg 


2 6 3-5 









D A 


Eg 


2 • 3 2 5 2 7 


D 4 






D 5 


Eg 


2 3 3 3 5 ■ 7 


A 3 






D 6 


Eg 


2 2 3 3 5 • 7 


AiAi 






D 7 


Eg 


2 3 3 3 5 









D s 


Eg 


3 3 5 









Eq 


Eg 


2 5 5 • 7 


A 2 






E 7 


Eg 


2 3 3-5 


Ai 






Eg 


Eg 


1 
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9 Eliminating root systems a priori 



There are 405844 root systems of rank n < 32 with no vectors of norm 1, 
corresponding to all direct sums of Ai,Di (i > A),Ee,Ej and Eg, where the 
order of the summands does not matter and the sum of the subscripts is at 
most 32. Since computing the number of embeddings r(i?i, Rj) can be time- 
consuming, we used the following congruences, due to Borcherds Q, to eliminate 
some root systems from consideration: 

Let roots(i?) denote the number of roots of R. If R is the root system of a 
32-dimcnsional even unimodular lattice then 

if R contains E%, then roots(i?) = (mod 24) 

if R contains E?, then roots(i?) = (mod 12) 

if R contains E 6 , then roots(-R) = (mod 6) 

if R contains D$, then roots(i?) = (mod 4) 

if R contains D?, then roots(i?) = (mod 8) 

if R contains Dg, then roots(i?) = (mod 8) 

if R contains D n ,n > 8, then roots(i?) = (mod 16) 

Also note that if a(R) = then m(R) = 0. If rank(i?) = 32 and det(i?) is 
not a perfect square then a(R) must be 0, so we eliminated those root systems 
as well. This left 135443 root systems. We ordered them so that dim(i?i) < 
dim(-Rj) if i < j, and so that det(-Rj) > det(i?j) if i < j and dim(i?i) = dim(i?j). 

As the matrix U with C/y = r(Ri,Rj) would still contain around 10 billion 
elements, we did not explicitly construct and invert it. Rather, we computed 
each element in the matrix when it was required for solving —Uv = w by 
back-substitution, with 

m{Ri) = r ^ R .j im- a(Ri) -^r(R j> R 4 )m(R j ) \ . 

If m(Rj) has already been computed to be 0, then the values r(Rj, R4) need 
not be computed since they make no contribution to this sum. 

Acknowledgements 

The author would like to thank Richard Borcherds for many helpful sugges- 
tions, Richard Fateman for advice on Lisp, and Xerox PARC for the use of its 
computers. This work was partially supported by grants from the NSF and the 
Royal Society. 



25 



References 

[1] R. Bacher, Unimodular lattices without nontrivial automorphisms, Inter- 
not. Math. Res. Notes 2 (1994) 91-95. 

[2] R. Bacher and B. B. Venkov, Reseaux entiers unimodulaires sans racines 
en dimension 27 et 28, Preprint No. 332, Inst. Fourier, Grenoble, 1996. 
Available at http://www-fourier.ujf-grenoble.fr/PREP/html/a332 

[3] E. Bannai, Positive definitive unimodular lattices with trivial automor- 
phism group, Mem. Amer. Math. Soc. 429 (1990) 1-70. 

[4] R. E. Borchcrds, The Leech lattice and other lattices, Ph.D. Dissertation, 
University of Cambridge, 1984. Available at arXiv:math. NT/9911195 
Much of this material also appears in [5] . 

[5] R. E. Borcherds, Classification of positive definite lattices, Duke Math. J. 
105 (2000), no. 3, 525-567. Available at arXiv: math. NT/99 12236 

[6] R. E. Borcherds, E. Freitag and R. Weissauer, A Siegel cusp form of degree 
12 and weight 12, J. Reine Angew. Math. 494 (1998) 141-153. 
Available at arXiv: math. AG/9805 132 

[7] J. H. Conway and N. J. A. Sloane, Low-dimensional lattices IV: the mass 
formula, Proc. R. Soc. Lond. A 419 (1988), 259-286. 

[8] J. H. Conway and N. J. A. Sloane, A note on optimal unimodular lattices, 
J. Number Theory 72 (1998), 357-362. 

[9] J. H. Conway and N. J. A. Sloane, Sphere Packings, Lattices and Groups, 
Springer- Verlag, NY, 3rd edition, 1998. 

[10] T. H. Cormen, C. E. Leiserson and R. L. Rivest, Introduction to Algorithms, 
MIT Press, Cambridge, MA, 1990. 

[11] N. D. Elkies, Lattices and codes with long shadows, Math. Res. Lett. 2 
(1995) 543-651. 

[12] N. D. Elkies, Mock-laminated lattices, preprint, 1996. 

[13] N. D. Elkies and B. H. Gross, The exceptional cone and the Leech lattice, 
Internat. Math. Res. Notices 14 (1996), 665-698. 

[14] M. R. Garey and D. S. Johnson, Computers and Intractibility: A Guide to 
the Theory of NP- Completeness, Freeman, SF, 1979. 

[15] H. Katsurada, An explicit formula for the Fourier coefficients of Siegel- 
Eisenstein series of degree 3, Nagoya Math J. 146 (1997), 199-223. 

[16] H. Katsurada, An explicit formula for Siegel series, Amer. J. Math. 121 
(1999), 415-452. 

[17] G. Kaufhold, Dirichletsche Reihe mit Funktionalglcichung in der Theorie 
der Modulfunktion 2: Grades, Math. Ann. 137 (1959), 454-476. 



26 



[18] M. Kervaire, Unimodular lattices with a complete root system, L'Enseign. 
Math. 40 (1994), 59-140. 

[19] O. D. King, Table of masses of even unimodular 32-dimcnsional lattices 
with any given root system. Available at arXiv: math. NT/0012231 

[20] Y. Kitaoka, A note on local densities of quadratic forms, Nagoya Math J. 
92 (1983), 145-152. 

[21] Y. Kitaoka, Local densities of quadratic forms and Fourier coefficients of 
Eisenstein series, Nagoya Math J. 103 (1986), 149-160. 

[22] Y. Kitaoka, Arithmetic of Quadratic Forms, Cambridge Tracts in Math., 
vol. 106, Cambridge Univ. Press, Cambridge, 1993. 

[23] Y. Kitaoka, Dirichlet series in the theory of quadratic forms, Nagoya Math 
J. 92 (1984), 73-84. 

[24] M. Kneser, Klasscnzahlen definitcr quadratischer Formcn, Arch. Math. 8 
(1957), 241-250. 

[25] H. Koch and B. B. Vcnkov, Uber ganzzahlige unimodulare euklidische Git- 
ter, J. Reine Angew. Math. 398 (1989), 144-168. 

[26] J. Leech, Notes on sphere packings, Canadian J. Math. 19 (1967), 251-267. 

[27] H. Maafi, Die Fourierkoeffizienten der Eisensteinrcihcn zwciten Grades, 
Mat. FysMedd. Dan. Vid. Selsk. 34 (1973), 1-13. 

[28] J. Milnor and D. Husemoller, Symmetric Bilinear Forms, Springer- Verlag, 
Berlin 1973. 

[29] G. Nebe and N. J. A. Sloane, A Catalogue of Lattices, published electroni- 
cally at http : //www. research, att . com/^nj as/lattices/ 

[30] H.-V. Nicmcicr, Definite quadratischc Formcn der Dimension 24 und 
Diskriminante 1, J. Number Theory 5 (1973), 144-178. 

[31] R. E. O'Connor and G. Pall, The construction of integral quadratic forms 
of determinant 1, Duke Math. J. 11 (1944), 319-331. 

[32] M. Peters, On even unimodular 32-dimensional lattices, Preprint SFB 
478, Mathematischen Instituts der Westfalischen Wilhelms-Universitat 
Minister, January 2001. 

[33] J. -P. Serre, A Course in Arithmetic, Springer- Verlag, NY 1973. 

[34] B. B. Venkov, The classification of integral even unimodular 24-dimcnsional 
quadratic forms. Trudy Matematicheskogo Instituta imeni V. A. Steklova 
148 (1978), 65-76. Also Chapter 18 of [9]. 

[35] G. L. Watson, Integral Quadratic Forms, Cambridge Univ. Press, Cam- 
bridge, 1960. 



27 



